Global existence and exponential growth for a viscoelastic wave 
equation with dynamic boundary conditions 
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Abstract 

The goal of this work is to study a model of the wave equation with dynamic boundary 
conditions and a viscoelastic term. First, applying the Faedo-Galerkin method combined with 
the fixed point theorem, we show the existence and uniqueness of a local in time solution. 
Second, we show that under some restrictions on the initial data, the solution continues to exist 
globally in time. On the other hand, if the interior source dominates the boundary damping, 
then the solution is unbounded and grows as an exponential function. In addition, in the absence 
of the strong damping, then the solution ceases to exist and blows up in finite time. 

Keywords: Damped viscoelastic wave equations, global solutions, exponential growth, blow up in 
finite time, dynamic boundary conditions. 
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1 Introduction 

We consider the following problem 



utt — Aw — aAut + 
u{x, t) = 0, 



g{t — s)Au{s)ds 



utt{x,t) 



du 



du 

u{x,0) = uo{x), 



{x,t) 



~ •5)7^1^) H — 7^ — (2;, t) + h [ut] 





ut{x,0) = ui{x) 



d. 



V 



X G 0, t > , 
X G To, t > 0, 
X G Ti, t > 0, 
X G , 



(1) 



where u = u{x,t) , t > , x € Q , A denotes the Laplacian operator with respect to the x variable, 
O is a regular and bounded domain of , {N > 1), dQ = Fg U Ti, mes(ro) > 0, Fq n Ti = 
and d/dv denotes the unit outer normal derivative, a is a positive constant, p > 2, h and g are 
functions whose properties will be discussed in the next section, uq , ui are given functions. 

Nowadays the wave equation with dynamic boundary conditions are used in a wide field of 
applications. See [24] for some applications. Problems similar to (1) arise (for example) in the 
modeling of longitudinal vibrations in a homogeneous bar in which there are viscous effects. The 
term Aut, indicates that the stress is proportional not only to the strain, but also to the strain 
rate, see [6] fore more details. 
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From the mathematical point of view, these problems do not neglect acceleration terms on the 
boundary. Such type of boundary conditions are usually called dynamic boundary conditions. 
They are not only important from the theoretical point of view but also arise in several physical 
applications. For instance in one space dimension and for g = 0, problem (1) can modelize the 
dynamic evolution of a viscoelastic rod that is fixed at one end and has a tip mass attached to its 
free end. The dynamic boundary conditions represents Newton's law for the attached mass, (see 
[5, 2, 8] for more details). In the two dimension space, as showed in [25] and in the references therein, 
these boundary conditions arise when we consider the transverse motion of a flexible membrane O 
whose boundary may be affected by the vibrations only in a region. Also some dynamic boundary 
conditions as in problem (1) appear when we assume that Q is an exterior domain of in which 
homogeneous fluid is at rest except for sound waves. Each point of the boundary is subjected to 
small normal displacements into the obstacle: this type of dynamic boundary conditions are known 
as acoustic boundary conditions, see [3] for more details. 

Littman and Markus [19] considered a system which describe an elastic beam, linked at its free 
end to a rigid body. The whole system is governed by the Euler-Bernoulli Partial Differential 
Equations with dynamic boundary conditions. They used the classical semigroup methods to 
establish existence and uniqueness results while the asymptotic stabilization of the structure is 
achieved by the use of feedback boundary damping. 

In [14] the author introduced the model 

utt - Uxx - utxx = 0, X e (0, L), t > 0, (2) 

which describes the damped longitudinal vibrations of a homogeneous flexible horizontal rod of 
length L when the end x = is rigidly fixed while the other end x = L is free to move with an 
attached load. Thus she considered Dirichlet boundary condition at x = and dynamic boundary 
conditions at x = L , namely 

uu{L,t) = -[u^ + ut:,]{L,t), t>0. (3) 

By rewriting the whole system within the framework of the abstract theories of the so-called B- 
evolution theory, the existence of a unique solution in the strong sense has been shown. An 
exponential decay result was also proved in [15] for a problem related to (2)-(3), which describe the 
weakly damped vibrations of an extensible beam. See [15] for more details. 
Subsequently, Zang and Hu [29], considered the problem 

Utt - P iux)xt - 1 = 0) X G (0, 1) , t > 

with 

u (0, t) = 0, p {u^\ + q K) (1, t) + kutt (1, t) = 0, i > 0. 

By using the Nakao inequality, and under appropriate conditions on p and g, they established both 
exponential and polynomial decay rates for the energy depending on the form of the terms p and 

q- 

Recently, the present authors have considered, in [11] and [12], problem (1) with g = and a 
nonlinear boundary damping of the form h (ut) = \ut\"^~'^ ut- A local existence result was obtained 
by combining the Faedo-Galerkin method with the contraction mapping theorem. Concerning the 
asymptotic behavior, the authors showed that the solution of such problem is unbounded and grows 
up exponentially when time goes to infinity provided that the initial data are large enough and 
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the damping term is nonlinear. The blow up result was shown when the damping is linear (i.e. 
m = 2). Also, we proved in [12] that under some restrictions on the exponents m and p, we can 
always find initial data for which the solution is global in time and decays exponentially to zero. 
These results had been recently generalized for a wide range of nonlinearities in the equation and 
in the boundary term: the authors proved the local existence and uniqueness by a sophisticated 
application of the non linear semigroup theory, see [13]. 

In the absence of the strong damping aAut and for Dirichlet boundary conditions on the whole 
boundary dQ, the question of blow up in finite time of problem (1) has been investigated by many 
authors. Messaoudi [23] showed that if the initial energy is negative and if the relaxation function 
g satisfies the following assumption 



then the solutions blow up in finite time. In fact this last condition has been assumed by other 
researchers. See for instance [16, 17, 21, 22, 26, 28]. 

The main goal of this paper is to prove the local existence and to study the asymptotic behavior 
of the solution of problem (1). 

One of the main questions is to show a blow-up result of the solution. This question is a difficult 
open problem, since in the presence of the strong damping term, i.e. when a ^ 0, the problem has 
a parabolic structure, which means that the solution gains more regularity. However, in this paper, 
we give a partial answer to this question and show that for a ^ and for large initial data, the 
solution is unbounded and grows exponentially as t goes to infinity. While for the case a = 0, the 
solution has been shown to blow up in finite time. 

The main contribution of this paper in this blow up result is the following: the exponential 
growth and blow-up results hold without making the assumption (4). In fact the only requirement 
is that the exponent p has to be large enough which is a condition much weaker than condition 
(4). Moreover, unlike in the works of Messaoudi and coworkers, we do not assume any polynomial 
structure on the damping term h{ut), to obtain an exponential growth of the solution or a blow up 
in finite time. 

This paper is organized as follows: firstly, applying the Faedo-Galerkin method combined with 
the fixed point theorem, we show, in Section 2, the existence and uniqueness of a local in time 
solution. Secondly, under the smallness assumption on the initial data, we show, in Section 3, 
that the solution continues to exist globally in time. On the other hand, in Section 4, we prove 
that under some restrictions on the initial data and if the interior source dominates the boundary 
damping then the L^-norm of the solution grows as an exponential function. Lastly, in Section 5, 
we investigate the case when a = and we prove that the solution ceases to exist and blows up in 
finite time. 

2 Preliminary and local existence 

In this section, we introduce some notations used throughout this paper. We also prove a local 
existence result of the solution of problem (1). 




(4) 



We denote 



H^^{n) = {ueH\n)/ur„=0}. 
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By (., .) we denote the scalar product in L^{Q) i.e. (n, v){t) = u{x, t)v{x, t)dx. Also we mean by 
||.||q the L'^{Q) norm for 1 < g < oo, and by ||.||5,ri the L'^{Ti) norm. 

Let T > be a real number and X a Banach space endowed with norm LP{0,T; X), 1 < 

p < oo denotes the space of functions / which are over (0, T) with values in X, which are 
measurable and \\f\\x G (0, T). This space is a Banach space endowed with the norm 

\\f\\LH0,T-,x) = llfrxdtY' . 

L°° {0,T; X) denotes the space of functions / : ]0,T[ ^ X which are measurable and \\f\\x £ 
(0,r). This space is a Banach space endowed with the norm: 



L°°{o,r;X) = ess sup ||/||x • 

0<t<T 

We recall that if X and Y are two Banach spaces such that X (continuous embedding), then 

LP (0, T; X) ^ LP{0,T;Y), l<p< oo. 
We will also use the embedding (see [1, Therorem 5.8]): 

2N 

H^^{n)^ LP{n), 2<p<p where p=\ if iV > 3, 



and also 

Hl^{^)^L'i{Ti),2<q<q where q 



+0O ifiV = l,2, 
+00 UN = 1,2. 



For 2 <m<q, let us denote V = H^^{n) D U''{Ti). 
We assume that the relaxation functions g is of class on M and satisfies: 

Vs G M, 5(s) > 0, and ^1 - ^ g (s) ds^ = I > . (5) 

Moreover, we suppose that: 

Vs > 0, ^'(s) < 0. (6) 
The hypotheses on the function h are the following: 

(HI) h is continuous and strongly monotone, i.e. for 2 < m < q, there exists a constant mo > 
such that 

{h{s) - h{v)) {s - v) > mo\s - vr , (7) 
(H2) there exist two positive constants Cm and Cm such that 

Cm\sr < Hs)s < Cm\sr, Vs G M . (8) 
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For a function u € 0, T], HY^{Q)j , let us introduce the following notation: 

{gou){t)= I g{t- s)\\Vu{s) -Vu{t)\\lds. 



Thus, when uec{[ 0, T], nC^ 0,r],L2(J7)^ such that ut G L'^ (o,T; H^^{n)^ , we have: 



d /•* 

- ((7 o ^x) (t) = g' {t - s) II Vn (s) - Vn 



Ids 



+-^[\\yu{t)\\l) I g{s)ds-2 I I g{t- s)Vu{s)Vut{t)dsdx 



n Jo 

{g ou){t)-2 / g{t- s)Vu (s) Vut (t) dsdx (9) 
Jo. Jo 

+|{l|Vu(t)||^ fg{s)ds\-g{t)\\Vuml. 



This last identity implies: 



j^j^g{t- s)Vu{s)Vut{t)dsdx = ]^[g' ^u){t) + ~^i^\Vu{t)\\l j^g{s)ds 



-\9{t)\\yu{t)\\l-\^^{gou){t) . (10) 



For ueC{[ f),T],Hl^{n)jr\C^[[ Q,T],L'^{n)j such that ut G L"^ {0,T; Hl^^{n)j , let us define the 
modified energy functional E by: 



= £;(t) = ^\\ut (t) \\l + ^\\ut (t) Wlr, + ^ (^l - 9{s)ds^ l|Vn(t)|| 



+^(50^)(t)-^||n(t)||^. (11) 

The following local existence result of the solution of problem (1) is closely related to the one 
we have proved for a slightly different problem in [11, Theorem 2.1], where no memory term was 
present. Let us sate it: 



Theorem 2.1. Assume that (5), (6) and (7) hold. Let 2 < p < q and max ^2, j < m < q. 

Then given uq G H^^{Q) and ui G L^(rj), there exists T > and a unique solution u of the problem 
(1) on (0, T) such that 



u G c([o,r],ifi,^(f])jnci([o,T],L2(o)j, 

Ut G L2(0,r;i?f„(17))nL™((0,T) xFi). 



Let us mention that Theorem 2.1 also holds for a = 0. The proof of Theorem 2.1 can be done 
along the same line as in [11, Theorem 2.1]. The main idea of the proof is based on the combination 
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between the Fadeo-Galerkin approximations and the contraction mapping theorem. However, for 
the convenience of the reader we give only the outhne of the proof here. 
For u G C{[0,T],H^^{n))nC^{[0,T],L^{n)) given, let us consider the following problem: 



X eCl, t> 0, 
X G To, t> 0, 
X eTi, t> 0, 
X G 0. 



(12) 



Vtt — — cxAvt + / g{t — s)Av{s)ds = \u\^ '^u, 
Jo 

v{x, t) = 0, 

Ov dv otdvt 

vtt{x,t) = - -^{x,t) - g{t - s) — {x,s)ds + ——{x,t) + h{vt) 
yov Jo ou 

v{x,Qi) = uo{x), vt{x,0) = ui{x) 

Definition 2.2. A function v{x,t) such that 

V G L~(0,r;/7f„(J7)) , 
Vt G L2(0,T;i7f„(f]))nL™ ((0,T)xri) , 
Vt G L'^{0,T;Hl^{Q))nL^{0,T;L\ri)) , 
Vtt G L°^{0,T;L\n))nL^{0,T;L\Ti)) , 
v{x,0) = tio(x), 
vt{x,0) = ui{x), 

is a generalized solution to the problem (12) if for any function io G H^^{^) n L"^{Ti) and ip G 
C^(0,r) with (p(T) = 0, we have the following identity: 



{\u\P-^u,w){t)^{t)dt 



T 



ivtt,w){t) + {Vv,Vw){t) 



+a{S/vt,Vw){t)ip{t)dt 



+ 



T 



git - s){Vv{s),Vw {t))ds 
vtt{t)wdr + h{vt) dT] dt. 



Lemma 2.3. Let 2 < p < q and 2 <m <q. Let uq e H'^{fl)nV, ui G -ff^(O), then for any T > 0, 
there exists a unique generalized solution (in the sense of Definition 2.2), v{t,x) of problem (12). 

The proof of Lemma 2.3 is essentially based on the Fadeo-Galerkin approximations combined 
with the compactness method and can be done along the same line as in [11, Lemma 2.2], we omit 
the details. 

In the following lemma we state a local existence result of problem (12). 



Lemma 2.4. Let 2 < p < q and max ^2, < m < q. Then given uq G H^^{i}) , ui G L^(r2) 

there exists T > and a unique solution v of the problem (12) on (0,T) such that 



V G c([o,r],ifi„(i7)jnci([o,T],L2(i7) 

Vt G L2(o,r;i^i^(0))nL-((0,r) xFi) 
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and satisfies the energy inequality: 



ht (t) M + ht it) lllr, + 1 - / 9{s)ds \\Vu (t) Wi + igou) (t) 



a / \\VvtiT)gdT + 




h{vt{cr, T))dadT 



< 




u{t)\^ '^u{T)vt{T)dTdx 



s Jn 



forO<s<t<T. 



Proof. We first approximate u e C([0, T],H^^ {Q)) D ([0, T],L'^{n)) endowed with the standard 

norm ||n|| = max ||nt(i)||2 + 11^(0 II /fMnii by a sequence {u'')k^N C C°°([0,T] x Q) hj a standard 

te[o,T] ^ ' 

convolution arguments (see [4]). Next, we approximate the initial data ui € L^(il) by a sequence 
{u'DkeN C C^{fi). Finally, since the space H^{n) nVn H^^^{n) is dense in H^^^in) for the H'^ 
norm, we approximate uq G H^^{Q) by a sequence (tio)fcGN C H^{Q.) nV H H^^{Q). 
We consider now the set of the following problems: 



Av^-aAv^+ [ g{t-s)Av''{s)ds 
Jo 



v^{x,t) = 0, 



^(a;,t)-y^ g{t - s) — {x,s)ds + -^{x,t) + 



advf 



h V 



^ v''{x,0) = ug, vf{x,0) 



u 



X eVL, t>0, 
X G To, t>0, 

X eTi, t > 0, 



(13) 

Since every hypothesis of Lemma 2.3 are verified, we can find a sequence of unique solution (ufc)fcgf^ 
of the problem (13). Our goal now is to show that {v'',vf)k£n is a Cauchy sequence in the space 

Yt ={ iv,vt)/vec{[o,T],H^^^in))ncmo,T],L^n)), 

vt G L2 (0, T; H^^ (f])) n ((0, T) x Ti) } 



endowed with the norm 



IIK^t)llYr = max ||?;dli + ^l|V^;||i +||i;t| 



0<t<T 



L™((0,T)xri 



0+ / w'^vtismds. 



(14) 



For this purpose, we set U = — u^' , V = — v^' . It is straightforward to see that V satisfies: 

( r* 

Vtt - Ay - aAVt + g{t- s)AV{s)ds = \u'' - l^-^n^' x G t > , 

Jo 

V{x,t) = X G Fo, t > 0, 



Vtt(x,t) 



-g^{x,t) - g{t - s) — {x,s)ds + + Kvt) - Hvt ) 



V{x, 0) = Uq - u^' , Vt{x, 0) = u'l 



,k' 



X G Fi, t>0, 
X G Jl. 
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We multiply the above differential equations by Vt, we integrate over (0, t) x Q, we use integration 
by parts and the identity (9) to obtain: 

I {\\Vt (t) Hi + \\Vt {t) ||2 + - / 5 W 11^^ (*)ll2) 

+aj^ WVVtWlds + ^ ^hiv1ix,T))-hivfix,T))) (v'[ix,T)-vfix,r))drdT 

= ^(ll^*(o)i + l|vy(o)||2 + ||y,(o)|||rJ 



Consequently, the above inequality together with (5), (6) and (7) gives 

\ {\\Vt {t) Hi + \\Vt (t) |||r, + I l|Vy ml) +aj'^ WVVtWlds + mo j'^ 11^*11™^'^^ 
< ^(ll^4(0)||i + ||VF(0)||i + ||y*(0)||i,rJ (16) 

Following the same method as in [11], we deduce that there exists C depending only on 17 and p 
such that: 

ll^lly, < C {\\Vtm\l + W^Vml + l|Vt(0)|||rJ + CT\\U\\y,. 

Since (no)fceN is a converging sequence in H^^ (fi), (u^)fcgN is a converging sequence in (17) 
and {u^)^^^ is a converging sequence in C ([0,r] ,iJpjj(Q)) n C-*^ ([0,T] ,L^(r2)) (so in also), we 
conclude that {v^ , is a Cauchy sequence in Yt- Thus (v^^v^) converges to a limit (u, ft) G Yr- 

Now by the same procedure used by Georgiev and Todorova in [10], we prove that this limit is a 
weak solution of the problem (12). This completes the proof of the Lemma 2.4. □ 

Proof of Theorem 2.1. In order to prove Theorem 2.1, we use the contraction mapping theorem. 
Indeed, for T > 0, let us define the convex closed subset of Yt: 

Xt = {{y,vt) G Yt such that t;(0) = tio,vt(0) = ui} . 

Let us denote: 

BuXXt) = {v e XT-MWr <R}^ 

the ball of radius R in Xt. Then, Lemma 2.4 implies that for any u G Xt, we may define v = ^ {u) 
the unique solution of (12) corresponding to u. Our goal now is to show that for a suitable T > 0, 
$ is a contractive map satisfying <I> {Br{Xt)) C Br{Xt). 
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Let u G Br{Xt) and v = ^ (u). Then for all t G [0, T] we have as in (15): 

t ^ 
\\vt\\2 + l\\^v\\l + \\vt\\l^^+2a j \\Vvt\\lds + c ll^^tCri^^^ 



(17) 



t 



< + llVuolli + |lni|||p^ + 2 j j \u{t)\'' '^u{T)vt{T) dxdr. 

n 

Using Holder's inequality, we can control the last term in the right hand side of the inequality (17) 
as follows: 

t t 

1 1 \u{t) r-'u{T)vt{T)dxdT< J ||n(r)||^-)(^_2)|bt(r) \\,^/^,^_^^^,^p_,^fT 
on 

2N 

Since p < — — -, we have: 

2N 2N 

< 



{3N - Np + 2{p-l)) - N -2 

Thus, by Young's and Sobolev's inequalities, we get for all 5 > there exists C{6) > 0, such that 
for all t G (0, T) 

t t 
j j \u (r) (r) vt (r) dxdr < C{5)tR^^P-^^ +5 j \\Vvt (r) Wldr. 
on 

Inserting the last estimate in the inequality (17) and choosing 5 small enough such that: 

1 



Thus, for T sufficiently small, we have llvHy^ < R. This shows that v G Br (Xt)- 

To prove that $ is a contraction, we have to follow the same steps (up to minor changes) as in 
[11]. We omit the details. Thus the proof of Theorem 2.1 is finished. □ 



Remark 2.5. Let us say that the hypothesis on m, max \ 2, -^^ ^ j < m < q, is made to pass to 
the limit in the nonlinear term, by the same way we have used in [11, Equation (2.28)]. 

3 Global existence 

In this section, we show that, under some restrictions on the initial data, the local solution of 
problem (1) can be continued in time and the lifespan of the solution will be [0, oo). 



Definition 3.1. Let 2 < p < q, max \^2, |+fr^j < m < q, uq e H^^{Q) and ui G L'^{Q). We 
denote by u the solution of (1). We define: 



= supjr > , -u = u{t) exists on [0, T] | 
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Since the solution u gYt (the solution is "regular enough"), from the definition of the norm given 
by (14) ! ^st us recall that ifT^ax < 00, then 

lim ||Vii(i)||2 + |lnt(t)|l2 = +00. 

t-'CTf-riax 

IfTmax < oc, we Say that the solution of (1) blows up and that Tmax is the blow up time. 
IfTmax = CO; we Say that the solution of (1) is global. 

In order to study the blow up phenomenon or the global existence of the solution of (1), for all 
< t < Tmax, we define: 



(18) 



I{t) = I{u{t)) = (^- j\{.s)d.^\\Vu{t)\\l + {g<>u){t)-\\u\\l, 

J{t) = J{u{t)) = ^(^l- l^g{s)dsy\Vn{t)\\l + ^{gou){t)-^\\u\\P. (19) 
Thus the energy functional defined in (11) can be rewritten as 

E{u{t)) = E{t) = J{t) + ]^\\ut\\l + \\\ut\\lr,. (20) 

As in [9, 27], we denote by B the best constant in the Poincare-Sobolev embedding H^^{Q,) ^ U'{Q.) 
defined by: 

B-^ = inf {||Vn||2 : u G HI^^{Q.), \\u\\p = l] . (21) 
For uq G H^oi^) , ui G -L^(il), we define: 

^(0) = ^Ikilli + ^lkilllri + ^l|V^xo||^-^||txo||^. 

The first goal is to prove that the above energy E (t) defined in (11) is a non-increasing function 
along the trajectories. More precisely, we have the following result: 



Lemma 3.2. Let 2 < p < q, max (^2, -^^^^ j < m < q, and u be the solution of (1). Then, for all 
t > 0, we have 



= ^{g'ou){t)-^g{t)\\Vu{t)\\l-a\\Vut\\l- 1^ h{ut)utdr 



< 



^{g'ou){t)-a\\Vut\\l- l_ hiut)utdr, VtG[0,T^ax). (22) 



Proof. Multiplying the first equation in (1) by ut, integrating over Q, using integration by parts we 
get: 

|{^IKIIi + ^lkll|r. + ^l|v-lli-^IK} 

- [ [ git- s)Vuis)\7ut{t)dsdx (23) 
Jn Jo 

= —a\\\7ut\\2 — / h{ut)utdT. 
A simple use of the identity (10) gives (22). This completes the proof of Lemma 3.2. □ 



10 



Lemma 3.3. Let 2 < p < q, max ^2, ^_^1_p ^ < m < q. Assume that (5) and (6) hold. Then given 
uq € H^^{Q,) , ui G L^(il) satisfying 



.£;(0)(p-2)/2^ < 1, 



I \l{p-2) ' ' J ' (24) 

t /K) >o, 

we have: 

i(u(t))>o, vtG[o,r^ax). 

Proof. Since / {uq) > 0, then by continuity, there exists T* < Tmax) such that 

/(t)>o, vtG[o,r*] 

which imphes that for all t € [0, T*], 

Jit) = h{t) + ^S^(^l- l\{s)dsy\Vu{t)\\l + {gou){t) 

- ^{[^- l^9{s)dsy\Vu{t)\\l + {gou){t)Y (25) 
By using (5), (6), (20) and (22), we easily get, for ah t G [0,r*] 

/||v^(i)||^ < -^Jit), 

p z 

< -^E{t)<^E{0). (26) 
p — 2 P — 2 

From the definition of the constant B in (21), we first get: 

WG[0,r*], \\u{t)\\l < BP\\Vu(t)r, . 

Since we have: 

yt G [o,r*] , B^vu{t)r, = _||Vn(t)||r' {l\\Vu{t)\\l) , 
by exploiting (26) and (24), we obtain, for all t G [0,r*]: 

||n(t)||^ < /3/(||Vn(t)i) 

< /^(l-^ 9is)ds^ ||Vtx(i)||i 

< (^1- l\{s)dsy\Vu{t)\\l 

Therefore, by using (18), we conclude that 

/(t)>0, VtG[0,T*]. 
Using the fact that E is decreasing along the trajectory, we get: 

BP ( 2v \(p-2)/2 

vo<t<r_,-^^^i^(t)j </3<i. 

By repeating this procedure, T* is extended to Tmax- D 
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Now, we are able to state the global existence theorem. 

Theorem 3.4. Let 2 < p < q, max ^2, < m < q. Assume that (5) and (6) hold. Then 

given uq € H^^{^1) , ui € L'^{0,) satisfying (24)- Then the solution of (1) is global and hounded. 
Proof. To prove Theorem 3.4, usmg the definition of Tmax, we have just to check that 

\\Vuit)\\l + \\utit)\\l 

is uniformly bounded in time. To achieve this, we use (19), (20), (22) and (26) to get 

EiO) > E{t) = J{t) + ^\\ut\\l + ^\\ut\\lr, 

> P^\\vn{t)\\l + ^\\utml. (27) 



Therefore, 

where C is a positive constant, which depends only on p. □ 



\Vu{t)\\l + \\ut{t)\\l<CE{0) 



4 Exponential growth for a > 

In this section we will prove that when the initial data are large enough, the energy of the solution 
of problem (1) defined by (11) grows exponentially and thus so the norm. 
In order to state and prove the exponential growth result, we introduce the following constants: 

B „-p/(p-2) ^ _/l n 2 ^ _fl 1\ 2 



Bi = j, a, = B-'''^-^\ Er = [^--^-yi, E,= \^--^-yi (28) 

Let us first mention that E2 < Ei. 

The following Lemma will play an essential role in the proof of the exponential growth result, and 
it is inspired by the work in [7] where the authors proved a similar lemma for the wave equation. 

First, we define the function 



(s) ds) \\Vu {t) \\l + {gou) (t) . (29) 



Let us rewrite the energy functional E defined by (11) as: 



E{t) = ^\\ut {t) i + (t) lllr^ + ^7 (*) " ^ l|n(t)||^ • (30) 



Lemma 4.1. Let 2 < p < q, max {^2, j <m<q. Let u be the solution of (1). Assume that 

E{0)<Ei and ||Vno||2 > ai- (31) 
Then there exists a constant 02 > ai such that 

(7(t))^/' >Q2, VtG[0,r^ax) (32) 

and 

\\u{t)\\p>Bia2, € [0,r„iax). (33) 
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Proof. We first note that, by (30), we have: 

E{t) > ijit)-l\\um; 

z p 

> hit)-^{i\\vum,r 

z p 

> 1^1 (t) - ^ {t)f' (34) 
2 p 

1 

= -a^ - = F (a) , 

I P 

1 /2 

where q = (7 (i)) ' .It is easy to verify that F is increasing for < a < ai, decreasing for a > oi, 
F (a) — > —00 as a ^ +00, and 

F{ai) = Ei, 

where ai is given in (28). Therefore, since E{0) < Ei, there exists 02 > ai such that F (02) = 
E{0). 

1 /2 

If we set ao = (7 (0)) , then by (34) we have: 

F{ao) <E{0) = F{a2), 

which implies that ao > a2- 

Now to establish (32), we suppose by contradiction that: 

(7(to))'/'<«2, 

for some to > and by the continuity of 7 (.) , we may choose to such that 

(7(to))'/'>ai. 

Using again (34) leads to: 

E {to) > F (7 (to)'/') > F{a2) = E{0). 

But this is impossible since for all t > 0, E{t) < E {0). Hence (32) is established. 
To prove (33), we use (30) to get: 

ll{t)<E{0) + -\\n{t)\\l . 
z p 



Consequently, using (32) leads to: 



But we have: 



^ll^(t)ll^ > ll{t)-E{0) 
> ^al-E{0). 



^al -E{0) = \x\ - F (02) = . 



Therefore (33) holds. This finishes the proof of Lemma 4.1. □ 
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The exponential growth result reads as follows: 
Theorem 4.2. Suppose that (5), (6) (8) hold. Assume that 

2 < m and max (m, 2/1) < p <p. 
Then, the solution of (1) satisfying 

E{Q) < E2 and \\Vuo\\^ > ai, (35) 
grows exponentially in the U' norm. 

Remark 4.3. It is obvious that for g = 0, we have Ei = E2, and Theorem 4-2 reduces to Theorem 
3.1 in [11]. 

Remark 4.4. In Theorem 4. 2, the condition 

{p/2) - 1 



g{s)ds < 



{p/2) - 1 + (l/2p) 

used in [16, 17, 21, 22, 23, 26, 28] is unnecessary and our result holds without it. 

Remark 4.5. Let us denote ci = — |^ — 2E2 {Bia2)~^ • Since we have seen that 02 > ai, using 
the definition of E2, we easily get ci > 0. This constant will play an important role in the proof of 
Theorem 4.2 

Proof of Theorem 4-2. We implement the so-called Georgiev-Todorova method (see [10, 22] and 
also [20]). So, we suppose that the solution exists for all time and we will prove an exponential 
growth. For this purpose, we set: 

J^{t) = E2-E (t) . (36) 

Of course by (31) and (22) and since E2 < Ei, we deduce that ^ is a non-decreasing function. 
So, by using (30) and, (36) we get successively: 

0<Jif{0)<J^{t)<E2-Eit)<Ei-l^{t) + -\\u{t)\\l. 

z p 

On one hand as F(ai) = Ei and V t > , 7(t) > 02 > a\, we obtain: 

£;i-^7(t) <F(ai)-^a? 

On the other hand, since 

F{ai) - -a^ = , 

2 p 

we obtain the following inequality: 

< Jif{0) <J^{t) <^\\u{t)\\P, Vt>0. (37) 



14 



For e small to be chosen later, and inspired by the ideas of the authors in [11], we then define the 
auxiliary function: 

^{t) = ^{t)+e I utudx + e I utudV + ^ \\V u\\l . (38) 

Let us remark that is a small perturbation of the energy. By taking the time derivative of (38), 
using problem (1), we obtain: 



a||Vnt||2 + y" h{ut)utdT + e\\ut\\\- e\\Vu\\\ 
+e + e llutllg Ti - ^ / h{ut)u{x,t)da 

+ Vu{t) g{t- s) Vii (s) dsdx. (39) 
Jn Jo 



By making use of (8) and the following Young's inequality 



XY < ^ + (40) 



X, y > 0, A > 0, ^, G IR+ such that 1/^ + Xjv = 1, then we get 



/ h{ut)udT < Cm \ut\ 



m—2 TT-i 

UtudL 



< C-m— II^Cr,+C™^^A-™/(™-i)||ntCr^. (41) 



Now, the term involving g on the right-hand side of (39) can be written as 



Vu{t,x) g{t- s)Vu{s,x)dsdx = \\Vu{t)\\U g (s) ds ] (42) 







+ / Vu{t,x) / g {t — s) {Vu {s, x) — Vu {t,x)) dsdx. 
Jn Jo 

On the other hand, by using Holder's and Young's inequalities, we infer that for all /i > 0, we get 

Vn(t, x) / g {t — s) {Vu {s,x) — Vu {t,x)) dsdx 
n Jo 

< [' g{t-s)\\Vu{t)\\2\\Vu{s)-Vu{t)\\2ds (43) 
Jo 

<^{gou){t) + ^(^j^g {s) d^ II Vu (t) 111. 

Inserting the estimates (41) and (42) into (39), taking into account the inequality (43) and making 
use of (8), we obtain by choosing /x = 1/2 and multiplying by I 



l^'it) > a/||V^zt||^ + /(c™-C™e^^A-"^/(™-i)) ||^xi||:;:_r,+e/||nt| 



+el \\u\\l + el \\ut\\l^^ - Cmsl— \\u\Q^^ (44) 



-|(ffon) it)-el\\Vuit)\\l 
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We want now to estimate the term involving ||w||^pj in (44). We proceed as in [11]. Then, we have 

\\u\ 

which holds for: 



im,ri ^ ^ \W\\H''{n) 



N N -1 

m>l and 0<s<l, s> > 0, 

2 m 

where C here and in the sequel denotes a generic positive constant which may change from line to 
line. 

Recalling the interpolation and Poincare's inequalities (see [18]) 



< Cll-ullp ||Vn||2, 



|Vn||2 , 
I ^ 

12 ■ 



we finally have the following inequality: 



|u||^,ri < C'll^llJ 'l|Vu||^. 



If s < 2/m, using again Young's inequality, we get: 



< C 



\u\ 



+ llVnl 



2\ 2 



(45) 



(46) 



for l/fM + 1/9 = 1. Here we choose 6 = 2/ms, to get /i = 2/(2 — ms). Therefore the previous 
inequality becomes: 



I ll'Ti ^ r~i 
l^llmA < C 



ip\ {2-ms)p 



+ llVull 



(47) 



Now, choosing s such that: 



we get: 



< s < 



2{p — m) 



m (p — 2) ' 

2m (1 - s) 

^ ^ < 1. 



(2 — ms) p 

Once the inequality (48) is satisfied, we use the classical algebraic inequality: 

z" <{z + l) < \l + ^{z + uj) , Vz>0, < z> < 1 , a;>0. 



(48) 



(49) 



to obtain the following estimate: 

( 



m(l — s)2 

p\T2^, < d(||n||^ + ^(0) 



< dnu\\l + J^{t)\ , Vt>0 



(50) 



where we have set d = 1 + 1/^(0). Inserting the estimate (50) into (46) we obtain the following 
important inequality: 



< c 



l^llp + ^l|Vu||^ + ^(i) 



(51) 
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Keeping in mind that 1 = 1 — g (s) ds, in order to control the term || Vu||2 in equation (44), we 
preferely use (as J^{t) > 0), the following estimate: 



I 11"^ ^ 



l'"llp + ^l|V^x||^ + 2^(t) 



which gives finally: 



\u\ 



m.Ti 



< c 



2^2 + 1 + - Ik 



V 



2 

2,ri 



+ /- 1 



Since 1 — g (s) ds > I, then we obtain from above 



2E2+[1 + 



m\p-\\ut 



2 

2,ri 



{g o u) (t) 



Now, inserting (53) into (44), then we infer that: 

l^'it) > al\\Vut\\l + l(^Cm-Cme'^^X-'^/^"'-^^ 



+e{l + ICr, 



-c 



m 



+eU- Cml—C 1 + 



m 

+e{Cml—C-- 
m 2 



2 + e [I + IQ 
2 



m. C \\ut\ 



2 

2,ri 



P 



igou) it) -el\\S/u(t) \\l - 2C„d — CE2 

m 



From (36), we get 



J^{t) < E2-\[l- j\{s)ds^ 



g{s)ds]\\Vu{t)\\l 



--{gou) {t) + -\\u 
1 p 



< E2-'-\\Vu{t)\\l-ligou)it) + -\\u 
lip 



This last inequality gives 



1 



I ||Vn(t)||2 >2[je{t)-E2 + -{gou) (t) 



u 



p 



Consequently, (54) takes the form 



\ m 



I II"* 

\MmT, 



m 



C]\\ut\\ 



2,ri 



+£ i^l + ICm—Cj WutWi + eil + lC„ 
2 A™ / 2 

+£0 Cml — Cil + - 

p m \ p 

+£ ( Cml—C - - + {gou){t) + 2£j^ (t) - 2Crael—CE2. 
m 2 I m 



Ml - 2eE2 



(52) 



(53) 



(54) 



(55) 



(56) 
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Now to estimate the terms involving ||n||^ and E2 in (56), we simply write: 



Then by using (33), we get: 

2 



l--]\\u\Z-2E2>c^\\u\Z, 



where ci > is defined in Remark 4.5. Thus, (56) becomes: 

l + lCm—C) \\ut\\l + e(l + lC„^—C) ||ut||2 

+e\ci- Cml—C 1 + - - ^C^el—CE2 (Biaa)"^ \ \\u\\l 
y m \ p J m J 

(\m ; \ / \m 

C^/— C -- + l]{gou){t) + 2e{.J^{t) + C„J—CE2 
ml) \ m 

Notice that since Z < 1, we first have : 

VA > 0, QJ — C-- + 1 > . 
m 2 

At this point, we pick A small enough such that: 

\m / 2 \ A™ 

ci - C„,l—C 1 + - - iCmel—CE2 (^102)^^ > . 
m \ p J m 

Once A is fixed, we may choose £ small enough such that 

c.„-C„e^A— /(™-i)>0, 
m 

^(0) > 0. 

Consequently, we end up with the estimate: 

^' (t) > m {htWl + htWlr, + \\u\\l + (t) + E2) , vt > o . 

Next, it is clear that, by Young's and Poincare's inequalities, we have: 



^{t)<^[j^{t) + \\ut\\l + \\ut\\l^^ + \\Vu\\l 
Since J^(t) > 0, then for all t > 0, we have: 



^||Vn||^< 11^11^ + £;2, 
z p 



for some 7 > 0. 
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Thus, the inequahty (59) becomes: 

^(t) <C J^{t) + \\ut\\l+\\ut\\lj^^ + \\u\\l + E2 for some C>0. (61) 
From the two inequahties (58) and (61), we finally obtain the differential inequality: 



dt 



> {t) for some f^i > 0. (62) 



An integration of the previous differential inequality (62) between and t gives the following 
estimate for the function ^: 

^ (t) > ^ (0) ef'K (63) 

On the other hand, from the definition of the function from inequality (37) and for small values 
of the parameter e, it follows that: 

^{t)<^\\u\\l. (64) 

From the two inequalities (63) and (64) we conclude the exponential growth of the solution in the 
L^-norm. □ 



5 Blow up in finite time for a = 

In this section, we prove that in the absence of the strong damping —Aut, (i.e. a = 0), the solution 
of problem (1) blows up in finite time that is it exists < T* < oo such that ||«(t)||p ^ oo as 
t T*. 

The blow up result reads as follows: 
Theorem 5.1. Suppose that (5), (6) and (8) hold. Assume that 

2<m and max{m,2/l) < p <p. 
Then, the solution of (1) satisfying 

E{0)<E2, |lVno||2>ai, (65) 
blows up in finite time. That is \\u {t)\\p — )■ oo as t — )• T* for some < T* < oo. 

Remark 5.2. The requirement m > 2 in Theorem 5.1 is technical but it seems necessary in our 
proof. The case m = 2 cannot be handled with the method we use here. But the same result can be 
shown for m = 2 by using the concavity method. See [12] for more details. 

Proof of Theorem 5.1. To prove Theorem 5.1, we suppose that the solution exists for all time and 
we reach to a contradiction. Following the idea introduced in [10] and developed in [20] and [27], 
we will define a function L which is a perturbation of the total energy of the system and which will 
satisfy the differential inequality 

^ > (t) , (66) 

where > 0. Inequality (66) leads to a blow up of the solution in finite time T* > L {0)~'^ 
provided that L (0) > 0. 



19 



To do so, we define the functional L as follows: 

L (t) = ^-"^(t) + e / utudx + e / UtudV, (67) 

where the functional Jif is defined in (36), a is satisfying 

^ • p-m p-2 m-2 \ 
\p(m— Ij 2p 2m J 

where a is defined later in (71) and e is a small positive constant to be chosen later. Taking the 
time derivative of L{t) and following the same steps as in the proof of Theorem 4.2, we get (instead 
of inequality (44)), for aU A > 0, 



IL' (t) > Ic^ (1 - a) (t) ll^iC^ri - C^e^^A-'"/^™-^) htWZr, + ht\ 



+el \\u\\l + el \\ut\\\^^ - Cmel— \\u\C^^^ (69) 
-^^{gou) {t)-d\\Vuit)\\l 

Next, for large positive M, we select A"'"/^"'"^) = M^'" (t). Then the estimate (69) takes the 
form: 

IL'it) > (^k„(l-a)-MC„,e^^) jr--(t)||nt||:;:ri+enhi|l2 

(i/f— (m— 1) 

+d Ml + d \\u,\\lr, - Cr^d^^H^^--^) \\u\\Z^^^ (70) 



-|(5on) {t)-d\\Vu{t)\\l 
Exploiting (37) and (45), we get: 

Thus, as in section 4, we have 



I sm 



|^||(l-.)m+ap(m-l) llVnll^'" < C 



\u\\lY ^ llVnIlp ^ 



ms6 



Choosing 9, and s exactly as in section 4 (with strict inequalities), we choose a that verifies: 

2-ms / 2m(l-s)\ „ 

a < — 1 - - — ^ = a. (71) 

- 2(m- 1) V {2-ms)pJ ^ ' 

The hypotheses on m and p ensure to have < o" < 1. 
Consequently, we get from above: 



wrx ' ^ +iivni 



(72) 
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Since, 



m(l-s) , 2 
^ '- + a (m - 1) I < 1, 



V 



2 — ms 



applying the algebraic inequality (49), we get: 



< d ||n||*; + jr(0) 



(73) 



Thus, (73) together with (72) leads to (see (53)): 



< Cd 



2^2 + 1 + - Ik 



2 - ll^tlllri -{gou) (t) 



(74) 



Inserting (74) into (70) and using (55), we obtain: 



IL'it) > {lcm{l-a)-MC^e^^^]j^-^{t)\\ut\i^^^ 



m 



I ]\/r—{m~l) \ , ^ 

+ e / [l + C^e Cd \ i^\\ut\\l + \\ut\\l^^j +2eje{t)-2eE2 



m 



(75) 



f 2 M-(""-i) / 2 

+ e ll Cml Cd 1 + - 

\ p m VP 



^j-(m.-l) 

\u\Z + Cmel Cd {g o u) [t] 



m 



j^-{m-l) / i' 

-2e C^l CdE2 + e[l--]{gou){t). 

m \ 2 



Writing again E2 = E2\\u\\p/\\u\\p and using again (33), we deduce that: 



IL' (t) > ( Icm. (1 - 0-) - MCme 



m — 1 



m 



m,ri 



e l\l + Cme Cd |j|iit||2 + ||nt||2^p^ 



m 



I + 2ejr (t) + 2C;„e/ CdE2 



m 

j^-(m-l) 



{2 M"(™~^) / 2 
/ 2E2 (^102)"^ - Cd 1 + - 1 - 4Cml 
p m \ p J m 

^-(m-l) 

+ e Cd(gou)(t). 

m 



CdE2 {Bia2y 
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Thus, using the definition of ci in Remark 4.5, we get: 

IL' (t) > (ic^ {I -a)- AfC„e^^^ jr^'^ (t) \\ut\C^^^ 



m 



+ e / 1 + Cme- 



+ e < ci - Cml 



j^-{m-l) 

m 

^-(m-l) 



Cd I \\ut\\l + |||p^ I + 2eJ^ (t) + 2C™e/ CdE2 



m 



m 



CdE2 {Bia2) 



J\,f-(m-l) 

+ e Cml Cd{gou){t) . 



m 



Since ci > 0, we choose M large enough such that: 



j^-{m-l) 



m 



-Cd{l + -] - ACml- 



^-(m-1) 



m 



CdE2 {Bia2)~^ > 0. 



Once M is fixed, we pick e smah enough such that 



Icm (1 - (t) - MCme > 

m 

and L (0) > 0. This leads to 

L' (t) > f? (llntll^ + ll-Uf III + (t) + ||n||^ + ^2 

for some r} > 0. 

On the other hand, it is clear from the definition (67), we have: 



J^{t)+ ( / ut udx + 



1 

l-cr 



{t) <C{e,cj) 

By the Cauchy-Schwarz inequality and Holder's inequality, we have: 

utudx < 



utudV 



1 

l-CT 



ufdx 



u^dx 



< C 



Ufdx 



|np dx 



(76) 



(77) 



where C is the positive constant which comes from the embedding (Q.) ^ (Q,). This inequality 
implies that there exists a positive constant Ci > such that: 



1 

l-CT 



Ufudx < Ci 



ItiP dx 



(1-<t)p 



ufdx 



2(l-a) 



Applying Young's inequality to the right hand-side of the preceding inequality, there exists a positive 
constant also denoted C > such that: 



1 

l-CT 



utudxj < C 



(1-ct)p 

|np dx ] + 



ufdx 



2(1-") 



(78) 
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for l/r + l/e = 1. We take 6* = 2(1 - a), hence r = 2 (1 - o") / (1 - 2o-), to get: 



1 

l-CT 



utudxj < C 



lup dx 



(l-2<7)v 



+ u^dx 
Jn 

Using the algebraic inequahty (49) with z = \\u\\^, d=l + 1/Jf{0), uj = Jf{0) and 
(the condition (68) on a ensuring that < < 1) we get: 

z" <d{z + jr(0)) <d{z + ,^{t)) . 

Therefore, there exists a positive constant denoted C2 such that for all t > 0, 



p{l-2a) 



1 

l-CT 



'^{t) + \\u{t)\\l+\\uty.jU2 



Ufudx I < C2 

Following the same method as above, we can show that there exists C3 > such that: 



(79) 



1 

l-CT 



UfudT j < C3 
Applying the inequality (51), we get: 



J^{t) + \\umZ,r, + \Wt it) 



\2Ti 



1 

l-CT 

UfudT ] < C4 



'j^it) + \\u{t)\\l + l\\Vu{t)\\l + \\ut (t)|||r^ 



Furthermore, inequality (60) leads to: 



UtudT 



V 



1 

— a 



<C5 



J^{t) + \\u{tW+\\ut{t) 



i2,ri 



+ E2 



Collecting (77), (79) and (80), we obtain: 

LT^ (t) < r/i {\\ut ml + \\ut\\l^^ + J^{t) + \\u{t)\\l + E2] , yt>0, 



(80) 



(81) 



for some 771 > 0. 

Combining (76) and (81), then, there exists a positive constant ^ > 0, as small as e, such that for 
all t > 0, 

L'{t)>iL^{t). (82) 

Thus, inequality (66) holds. Therefore, L{t) blows up in a finite time T*. 

On the other hand, from the definition of the function L{t) and using inequality (37), for small 
values of the parameter e, it follows that: 

l-a 

(83) 



L{t)<K[\\u{t)\\l 

where k is a positive constant. Consequently, from the inequality (83) we conclude that the norm 
||?x(i)||p of the solution n, blows up in the finite time T*, which implies the desired result. This 
completes the proof of Theorem 5.1. □ 
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